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Abstract. This paper is an introduction to some of the main present issues in the 
theory of structural glasses. After recalling a few experimental facts, it gives a short 
account of the analogy between fragile glasses and the mean field discontinuous spin 
glasses. The many valley picture is presented, and a brief account of recent attempts 
to obtain quantitative results from first principle computations is summarised. 



1. Introduction 

When quenched fast enough so that it avoids the crystallisation transition, almost any 
liquid becomes a glass This means that the density profile is not flat as in a liquid, 
it contains some peaks as in a crystal, but these peaks are not located on the nodes of 
a periodic or quasi periodic lattice. The understanding of such amorphous 'solid' states 
has been recognised for a long time as a major question in condensed matter physics. 
The sentence by Phil Anderson: "... there are still fascinating questions of principle 
about glasses and other amorphous phases..." ||, written nearly thirty years ago, was 
once again visionary in that it foresaw the wonderful developments on glassy systems, 
and particularly on spin glasses. The progress has been particularly difficult in these 
area, and in particular as far as structural glasses are concerned. 

2. Mathematics 

The first question which comes to mind is whether the glass is a new state of matter. It 
is not distinguished by any obvious symmetry (a not-obvious symmetry will be discussed 
later) from the liquid state, and one might think (as many people do) that the density 
profile would actually become flat on time scales longer than the experimental ones: the 
glass would just be a liquid with a long relaxation time. 

From the statistical physics point of view, one wants to start from a microscopic 
Hamiltonian. The simplest situation is that of iV point-like particles in a volume V, 
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with a pair interaction potential 



(1) 



i<j 



A simple case is that of homogeneous systems where Vtj is independent of i and j, and can 
be for instance either a hard sphere potential, a 'soft sphere' potential (Vij(r) = A/r 12 ), 
or a Lennard- Jones potential (V^-(r) = A/r 12 — B/r 6 ). Also much studied numerically 0, 
because the crystallisation is more easily avoided, are the binary mixtures where there 
are two types of particles: each particle i has G {±1} and V^-(r) = V^ i£ . (r), where 

V++, V , and V + - = V- + are three potentials of the same type as before, but with 

different A, B parameters corresponding to particles + and — having different radii. 

Does there exist, in any such case, an independent state of matter which is the 
glass state? Does it exist as a long-lived metastable state (like the diamond phase of 
carbon)? Nobody knows the rigorous mathematical answer to these questions. Actually 
much simpler related questions are unanswered (e.g. proving the existence of a spin glass 
phase in a finite dimensional short range system pj), or have taken many efforts to solve 
(e.g. proving Kepler's conjecture that the densest three dimensional packing of hard 
spheres is the fcc/hcp lattice ||). 

3. Experiments 

Experimentally, the liquid falls out of equilibrium on experimental time scales, and 
becomes a 'glass', at a temperature T g called the glass temperature Jl|]. This glass 
temperature is conventionnally defined as the one at which the relaxation time r of 
the liquid, as obtained e.g. from viscosity or from susceptibility measurements, becomes 
of the order of 10 3 seconds. Angell's plot of log (t/1s) versus T g /T allows to distinguish 
several types of behaviour (fig. [1]). So called strong glasses like Si0 2 have a typical 
Arrhenius behaviour with one well defined free energy barrier. On the other hand, some 
glasses, called fragile, show a dramatic increase of the relaxation time when decreasing 
temperature which is much faster than Arrhenius: the typical free energy barrier thus 
increases when T decreases. This implies a collective behaviour involving more and 
more particles. An increase of the dynamical correlation, characteristic of the mobile 
particles (rather than the more natural correlation of frozen particles), has been found 
in recent simulations [|7| . A popular fit of the relaxation time versus temperature is the 
Vogel Fulcher one, 



which would predict a phase transition at a temperature T VF which is not accessible 
experimentally (while staying at equilibrium). The more fragile the glass, the closer is 
Typ to T g , while strong glass formers have a Typ close to zero. 

Another interesting experimental signature is that of the specific heat. When one 
cools the liquid slowly, at a cooling rate T = —dT/dt, it freezes into a glass at a 
temperature which decreases slightly when F decreases. When this freezing occurs, 




(2) 




Figure 1. The left-hand figure shows the behaviour of the logarithm ( in base 10) of 
the relaxation time in seconds, versus T g /T, for two extreme cases of glass formers. The 
dashed line is from Ge02, a strong glass former which has an Arrhenius like behaviour; 
The full line is the measurement from OTP, which is a fragile glass former with a 
relaxation time diverging much faster than an Arrhenius law. The right-hand figure 
shows the configurational entropy S c of OTP (in J K~ 1 mol~ 1 ) versus temperature 
(in K). The configurational entropy, defined as the difference between the entropy 
of the supercooled liquid and that of the crystal, is measured through an integral 
of the specific heat difference. The squares are the experimental values. The glass 
temperature is T g = 246K. The full line is a fit to the equilibrated data, of the 
type S c — Sao(l — Tk/T). The Kauzmann temperature is Tjc = 204K, the fusion 
temperature is 331-ftT. The data are taken from Q. 



the specific heat jumps downward, from its value in the equilibrated supercooled liquid 
state to a glass value which is close to that of the crystal. From the specific heat, 
one can compute the entropy. The configurational entropy, defined experimentally as 
the difference S c = Su q — S crys tai, behaves smoothly in the supercooled liquid phase, 
until the system becomes a glass (see figjl]). It was noted by Kauzmann long ago 
that, if extrapolated, S C (T) vanishes at a finite temperature Tr. If cooled more slowly, 
the system follows the smooth S C (T) curve down to slightly lower temperatures, but 
then freezes again. One can wonder what could happen at infinitely slow cooling. As 
a negative S c does not make sense (except for pure hard spheres, where there is no 
energy), something must happen at temperatures above T K . The curve S C (T) could 
flatten down smoothly, or there might be a phase transition, which in the simplest 
scenario would lead to S C (T) = at T < T#. This idea of an underlying "ideal" 
phase transition, which could be obtained only at infinitely slow cooling, receives some 
support from the following observation: the two temperatures where the extrapolated 
experimental behaviour has a singularity, T VF and T K , turn out to be amazingly close to 
each other (see the table below) . The first phenomenological attempts to explain this 
fact originate in the work of Kauzmann || , and developed among others by Adam, Gibbs 
and Di-Marzio [[J, which identifies the glass transition as a 'bona fide' thermodynamic 
transition blurred by some dynamical effects. 
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If there exists a true thermodynamic glass transition at T = T# = Tvf-, it is a 
transition of a strange type. It is of second order because the entropy and internal 
energy are continuous. On the other hand the order parameter is discontinuous at the 
transition, as in first order transitions: the modulation of the microscopic density profile 
in the glass does not appear continuously from the flat profile of the liquid. As soon as 
the system freezes, there is a finite jump in this modulation (A more precide definition 
of the order parameter will be given below). 



Comparison of T# and Typ in various glass-formers (from 0) 



Substance 


Tk(K) 


Tvf(K) 


T g (K) 


o-terphenyl 


204.2 


202.4 


246 


salol 


175.2 




220 


2-MTHF 


69.3 


69.6 


91 


n-propanol 


72.2 


70.2 


97 


3-bromopentane 


82.5 


82.9 


108 



4. A mean field spin glass analogy 

A totally different class of systems where such a lst-2nd order type transition was found, 
and studied in great details, is a certain category of mean field spin glasses. A few years 
after the replica symmetry breaking (RSB) solution of the mean field theory of spin 
glasses PH[ , it was realized that there exists another category of mean-field spin glasses 



where the static phase transition exists and is due to an entropy crisis [0. These are 
now called discontinuous spin glasses because their phase transition, although it is of 



second order in the Ehrenfest sense, has a discontinuous order parameter [T2j . Another 
name often found in the literature is 'one step RSB' spin glasses, because of the special 
pattern of symmetry breaking involved in their solution. These are spin glasses with 
infinite range interactions involving a coupling between triplets (or higher order groups) 
of spins. The simplest among them is the random energy model, which is the p — ► oo 
limit version of the p-spin models described by the Hamiltonian 

H = — ^ ^ii...i p Sj 1 ...Sj p (3) 

iK...<i p 

where the J's are (appropriately scaled) quenched random couplings, and the spins can 
be either of Ising or spherical type [1^, [13], [L4| . 

The analogy between the phase transition of discontinuous spin glasses and the 
thermodynamic glass transition was first noticed by Kirkpatrick, Thirumalai and 
Wolynes in a series of inspired papers of the mid-eighties [[13] . While some of the basic 



ideas of the present development were around at that time, there still missed a few 
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C(t) 




Figure 2. The left-hand figure shows the schematic behaviour of the correlation 
function found in mean field discontinuous spin glasses and observed in structural 
glasses. The typical two-step relaxation consists of a fast (3 relaxation leading to a 
plateau, followed by a a relaxation from the plateau, whose typical time scale increases 
rapidly when T decreases, and diverges at T = T c which is equal to the mode-coupling 
transition temperature. The right-hand figure shows the behaviour of the relaxation 
time versus temperature. The right-hand curve is the prediction of mode-coupling 
theory without any activated processes: it is a mean field prediction, which is exact 
for instance in the discontinuous mean-field spin glasses |l8j . The left-hand curve is the 
observed relaxation time in a glass. The mode coupling theory provides a quantitative 
prediction for the increase of the relaxation time when decreasing temperature, at 
high enough temperature (well above the mode coupling transition T c )jl^, ^(J. The 
departure from the mean field prediction at lower temperatures is usually attributed 
to 'hopping' or 'activated' processes, in which the system is trapped for a long time in 
some valleys, but can eventually jump out of it. The ideal glass transition, which takes 
place at Tr-, cannot be observed directly since the system becomes out of equilibrium 
on laboratory time scales at the 'glass temperature' T g . 

crucial ingredients. On one hand one needed to get more confidence that this analogy 
was not just fortuitous. The big obstacle was the existence (in spin glasses) versus the 
absence (in structural glasses) of quenched disorder. The discovery of discontinuous spin 
glasses without any quenched disorder |15], [16], [T7j provided an important new piece of 
information: contrarily to what had been believed for long, quenched disorder is not 
necessary for the existence of a spin glass phase (but frustration is). 

It is important to analyse critically this analogy from the point of view of the 
dynamical behaviour. In discontinuous mean field spin glasses there exist a dynamical 
transition temperature at a temperature T c which is larger than the equilibrium 
transition T^- When T decreases and gets near to T c , the correlation function relaxes 
with a characteristic two step forms: a fast (3 relaxation leading to a plateau takes place 
on a characteristic time which does not grow, while the a relaxation from the plateau 
takes place on a time scale which diverges when T —>■ T c (see fig. |2|). This dynamic 
transition is exactly described by the schematic mode coupling equations. 
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However the existence of a dynamic relaxation at a temperature above the true 



thermodynamic one is possible only in mean field, and the conjecture |13[ is that in a 
realistic system like a glass, the region between Tk and T c will have instead a finite, but 
very rapidly increasing, relaxation time, as explained in fig. 0. A similar behaviour has 
been found in finite -size mean field models |2l| 



Another very interesting dynamical regime is the one where the system is out of 
equilibrium (T < T g ). Then the system is no longer stationnary: it ages. This is well 
known for instance from studies in polymeric glasses. If one measures the response of 
your favorite plastic ruler to some stress, it will behave differently depending on its 
age. Schematically, new relaxation processes come into play on a time scale comparable 
to the age of the system: the older the system, the longer the time needed for this 
"aging" relaxation to take place. Recent years have seen important developments 



on the out of equilibrium dynamics of the glassy phases p^| , initiated by the exact 
solution of the dynamics in a discontinuous spin glass by Cugliandolo and Kurchan 
It has become clear that, in realistic systems with short range interactions, the 



pattern of replica symmetry breaking can be deduced from the measurements of the 



violation of the fluctuation dissipation theorem ||24j| . These measurements are difficult. 
However, numerical simulations performed on different types of glass forming systems 
have provided an independent and spectacular confirmation of their 'one step rsb' 
structure [[25|, [26|, ^7|, on the (short) time scales which are accessible. Experimental 
results have not yet settled the issue, but the first measurements of effective temperatures 
in the fluctuation dissipation relation have been made recently |[29|| . 

To summarize, the analogy between the phenomenology of fragile glass formers and 
discontinuous mean field spin glasses accounts for: 

• The discontinuity of the order parameter 

• The continuity of the energy and the entropy 

• The jump in specific heat (and the sign of the jump) 

• Kauzmann's "entropy crisis" 

• The two steps relaxation of the dynamics and the succes of Mode Coupling Theory 
at relatively high temperatures. 

• The aging phenomenon and the pattern of modification of the fluctuation 
dissipation relation in the low temperature phase 

5. A lesson from mean field: many valleys 

The successes of the above analogy suggest to have a closer look at the mean field models 
in order to understand, at least at the mean field level, what are the basic ingredients at 
work in the glass transition. In mean field spin glasses, at temperatures T K < T < T c , 
the phase space breaks up into ergodic components which are well separated, so-called 



free energy valleys or TAP states fl30| , plf] . Each valley a has a free energy F a and a free 
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Sc(f) 




Figure 3. Qualitative shape of the complexity versus versus free energy in mean field 
discontinuous spin glasses. The whole curve depends on the temperature. The saddle 
point which dominates the partition function, for m constrained replicas, is the point 
/* such that the slope of the curve equals m/T (for the usual unreplicated system, 
m = 1). If the temperature is small enough the saddle point sticks to the minimum 
/ = fmin and the system is in its glass phase. For m = 1, this equilibrium phase 
transition happens at T = Tk- 



energy density f a = F a /N. The number of free energy minima with free energy density 
/ is found to be exponentially large: 

Af(f,T,N)nexv(NX(f,T)), (4) 

where the function E is called the complexity. The total free energy of the system, $, 
can be well approximated by: 



e 



P^ Ee W) = / d/ exp(iV[E(/,T) -/?/]) , (5) 



ma ;r 



where f3 = 1/T. The minima which dominate the sum are those with a free energy 
density /* which minimizes the quantity $(/) = f — T£(/, T). At large enough 
temperatures the saddle point is at / > f m i n (T). When one decreases T the saddle 
point free energy decreases (see fig.|3|, with m — 1). The Kauzman temperature T# is 
that below which the saddle point sticks to the minimum: /* = f m i n (T). It is a genuine 



phase transition |TT|, [12|, |13fl . However because T c > T K , the phase space is actually 
separated into non ergodic components (valleys) at T < T c (actually there exist some 
non ergodic components also above T c , but they are not felt by the system when starting 
from random initial conditions |[32| ). The total equilibrium free energy is analytic at T c : 
in spite of the ergodicity breaking, the system has the same free energy as that of the 
liquid, as if transitions were allowed between valleys. 

What remains of this mean field picture in finite dimensional glasses? When one 
decreases the temperature, there is a well defined separation of time scales between 
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the a and the (3 relaxations, which suggest to consider the dynamical evolution of 
system in phase space as a superposition of two processes: an intravalley (relatively 
fast) relaxation, and an intervalley (slow, ang getting rapidly much slower when one 
cools the system) hopping process. 

One popular way of making this statement more precise, allowing to study it 
numerically, is through the use of inherent structures (IS) ||33|| . Given a configuration of 
the system, characterized by its phase space position x = {xi, ...,x*n}, the corresponding 
inherent structure s(x) is another point of phase space which is the local minimum of 
the Hamiltonian which is reached from this configuration through a steepest descent 
dynamics. IS are easily identified numerically. A given trajectory x(t) of the glass 
through phase space maps onto the corresponding trajectory s(t) in the space of inherent 



structures. Looking at the dynamical evolution in the space of IS |34| makes the valley 
structure slightly more apparent, since one gets rid of the small thermal excitations 
around each valley minimum. Calling T> s the set of those configurations which are 
mapped to the coherent structure s, a natural definition of the IS entropy density, £j s , 
is NH is (T) = — J2 S P( s ) ^(/^(s)), where the weight of the inherent structure s is 

P{s) = Z{s)IY d Z{b) ; Z(s)= [ dxexp(-(3H(x)) . (6) 

In a system with short range interactions, it is reasonable to expect that one may 
have two distinct IS which differ by a local rearrangement of a finite number of atoms. 
It is then easy to show that the slope of configurational entropy versus free energy is 



infinite around j m i n p5 |, which does not agree with the general scenario, except if the 
Kauzmann temperature vanishes. This problem is due to te fact that IS are too simple 
objects, which cannot be identified with the free energy valleys. The difference is very 



easily seen in spin systems ||36|| : IS are nothing but configurations which are stable to 
one spin flip. Zero temperature free energy valleys, defined as TAP states, are stable 
to the flip of any arbitrarily large number k of spins (but the limit iV — > oo must be 
taken before the limit k — > oo). In continuous systems, the generalization is clear: IS 
are local minima of the energy, so that any infinitesimal move of the positions of all N 
particles raises the energy. Let us generalize the notion of a minimum as follows: define 
a k-th order local minimum as a configuration of particles such that any infinitesimal 
move of all N particles, together with a move of arbitrary size of k particles, raises 
the energy. The limit k — > oo gives the proper definition of a zero temperature free 



energy valley. The proper definition at finite temperature is slightly more involved B7 



Let us summarize it here briefly. Given two configurations x and y we define their 
overlap as before as q(x,y) = — l/NJ2ik=i n w(xi ~ Uk), where w(x) = —1 for x small, 
w(x) = for x larger than the typical interatomic distance. We add an extra term to 
the Hamiltonian: we define 

exp(-NpF(y,e)) = J dx exp(-if (x) + (3eNq(x, y)), 

F(e) = (F(y,e)}, ' (7) 
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where (f(y)) denotes the average value of / over equilibrium configurations y 
thermalized at temperature (3~ l . Taking the thermodynamic limit before the limit e — ► 
allows to identify the valley around any generic equilibrium configuration y J37j p9[] . 

In a nutshell, two configurations which differ by the (arbitrarily large) displacement 
of a finite number of atoms are in the same thermodynamic valley. This definition of the 
valleys also suffers from some difficulties: Nucleation arguments then forbid the existence 
of a non-trivial complexity versus free energy curve in a finite dimensional system. The 
solution consists in noticing that there exist many more metastable valleys, which have 
a finite but very long lifetime. These can be identified by taking the Legendre transform 
W(q) of the free energy F(e): 

-OF 

W(q)=F(e) + eq ; <7 = • (8) 



Analytic computation in mean field models [37], as well as in glass forming liquids using 



the replicated HNC approximation ]35|, show that W(q) is minimal at q — 0, but has 
a secondary minimum at a certain q = qEA-, in the temperature range Tk < T < T c . 
The behaviour around this second, metastable, minimum corresponds to phenomena 
that can be observed on time scales shorter than the lifetime of the metastable state. 
The thermodynamic configurational entropy is the value of the potential W(q) at the 
secondary minimum with q ^ f3"7fl , and it can be defined only if the minimum does 
exist (i.e. for T < T c ). Of course the secondary minimum for T > Tjt is always in the 
metastable region. However if one would start from a large value of e and would decrease 
e to zero not too slowly, the system would not escape from the metastable region and 
one obtains a proper definition of the thermodynamic configurational entropy in this 
region T > Tk- In a similar way one could compute q(e) in the region (e > e c ) where the 
high q phase is thermodynamically stable and extrapolate it to e — > 0. The ambiguity 
in the definition of the thermodynamic configurational entropy at temperatures above 
Tk becomes larger and larger when the temperature increases. It cannot be defined for 
T>T C . 

6. Beyond the analogy: first principles computation 

In recent years, it has become possible to go beyond the simple analogy between 
structural glasses and mean field discontinuous spin glasses. One can actually use the 
concepts and the techniques which are suggested by this analogy in order to start a 
systematic first principles study of the glass phase [fy], [|2| |39| . So far we have focused 



onto the equilibrium study of the low temperature phase. One main reason is that 
the direct study of out of equilibrium dynamics is more difficult, and that one might 



be able to make progress by a careful analysis of the landscape ||40|| . The strategy is 
to assume that there exists a phase transition, and that it is of the same type as the 
one found in discontinuous mean field spin glasses. Within this framework, one tries 
to compute the properties of the glass phase. This involves several quantities like the 
Kauzmann temperature, the radius of the cage which confine the particles in the glass 
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phase, the configurational entropy etc... The validity of the scenario is checked from 
the comparison of various predictions with numerical simulations of well equilibrated 
systems. 

The first task is to define an order parameter. This is not trivial in an equilibrium 
theory where we have no notion of time persistent correlations. The best way is to 
introduce two copies of the system, with a weak interaction. The two sets of particles 
have positions Xi and yt respectively, the total Hamiltonian is 

E = E ( v ( x i- x j) + v (yi-yj)) + e J2 w ( x i-yj) ( 9 ) 

where we have introduced a small attractive potential w(r) between the two systems. 
The precise shape of w is irrelevant, insofar as we shall be interested in the limit e — > 0, 
but its range should be of order or smaller than the typical interparticle distance. The 
order parameter is then the correlation function between the two systems: 

9x V {r) = Jim lirn^ ^ < 5 ( x i ~ Vi ~ r ) > ( 10 ) 

In the liquid phase this correlation function is identically equal to one, while it has 
a nontrivial structure in the glass phase, reminiscent of the pair correlation of a 
dense liquid, but with an extra peak around r ~ 0. Let us notice that we expect a 
discontinuous jump of this order parameter at the transition, in spite of its being second 
order in the thermodynamic sense. The existence of a non trivial order parameter is 
associated with the spontaneous breaking of a symmetry: For e = 0, with periodic 
boundary conditions, the system is symmetric under a global translation of the x 
particles with respect to the y particles. This symmetry is spontaneously broken in 
the low temperature phase, where the particles of each subsystem tend to sit in front of 
each other. 

Generalizing this approach to a system of m coupled replicas, sometimes named 
'clones' in this context (the order parameter used only m = 2), provides a wonderful 
method for studying analytically the thermodynamics of the glass phase [^, In the 



glass phase, the attraction will force all m systems to fall into the same glass state, so 
that the partition function is: 

Z m = Y, e-P Nmf «W = f fmax df exp (jvTOj T) _ m(3f]) (n) 



In the limit where m — > 1 the corresponding partition function Z m is dominated by the 
correct saddle point /* for T > Tk- The interesting regime is when the temperature is 
T < Tk, and the number m is allowed to become smaller than one. The saddle point 
f*(m,T) in the expression ( |TT| ) is the solution of dT>(f,T)/df = m/T. Because of the 
convexity of S as function of /, the saddle point is at / > f m i n (T) when m is small 
enough, and it sticks at /* = f m in(T) when m becomes larger than a certain value 
m = m*(T), a value which is smaller than one when T < Tk (see fig. |||). The free 
energy in the glass phase, F(m = 1,T), is equal to F(m*(T),T). As the free energy is 
continuous along the transition line m = m*(T), one can compute F(m*(T),T) from 
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the region m < m*(T), which is a region where the replicated system is in the liquid 
phase. This is the clue to the explicit computation of the free energy in the glass 
phase. It may sound a bit strange because one is tempted to think of m as an integer 
number. However the computation is much clearer if one sees m as a real parameter 
in flllf ). As one considers low temperatures T < Tk the m coupled replicas fall into 
the same glass state and thus they build some molecules of m atoms, each molecule 
being built from one atom of each 'colour'. Now the interaction strength of one such 
molecule with another one is basically rescaled by a factor m (this statement becomes 
exact in the limit of zero temperature where the molecules become point like). If m is 
small enough this interaction is small and the system of molecules is liquid. When m 
increases, the molecular fluid freezes into a glass state at the value m = m*(T). So our 
method requires to estimate the replicated free energy, F(m, T) = —\og(Z m )/(PmN), in 
a molecular liquid phase, where the molecules consist of m atoms and m is smaller than 
one. For T < Tk, F(m, T) is maximum at the value of m = m* smaller than one, while 
for T > Tk the maximum is reached at a value m* is larger than one. The knowledge 
of F m as a function of m allows to reconstruct the configurational entropy function 
Sc(f) at a given temperature T through a Legendre transform, using the parametric 



representation (easily deduced from a saddle point evaluation of ([11])) HI 



d[mF(m,T)} _ m 2 dF(m,T) 

am 1 am 

The Kauzmann temperature ('ideal glass temperature') is the one such that 
i7i*(Tk) = 1. For T < T# the equilibrium configurational entropy vanishes. Above 
Tk one obtains the equilibrium configurational entropy E(T) by solving (]i~2|) at m — 1. 

This gives the main idea which allows to compute the free energy in the glass 
phase, at a temperature T < Tk, from first principles: it is equal to the free energy of 
a molecular liquid at the same temperature, where each molecule is built of m atoms, 
and an appropriate analytic continuation torn — m*(T) < 1 has been taken. The whole 
problem is reduced to a computation in a liquid. This is not trivial, and requires to 
develop some specific approximations. I shall not elaborate on that here, but refer the 



reader to the original papers |42|, [45], fig, |47]| . The basic idea of the approximation is 
that the size of the molecules is directly related to the thermal wandering of an atom 
in its cage. Therefore at low temperatures one can use some small cage approximation, 
it is natural to write the partition function in terms of the center of mass and relative 
coordinates {r iy uf\ y with = + it" and J2 a u i = 0> an d to expand the interaction in 
powers of the relative displacements u. Keeping only the term quadratic in u (harmonic 
vibrations of the molecules), and integrating over these vibration modes, one gets the 
"harmonic resummation" approximation where the partition function is given by: 

Z m = Zl J dr exp (-(3mH(r) - ^^^Tr log m\ (13) 

-JVd(m-l) 



where Z^ = m Nd / 2 y/27rT /N\, and the matrix M, of dimension dN x dN, is given 
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by: 

d 2 H(r) 

M mw = * r ^ r u = s a E v A r i - r k) - v^in - r 3 ) (14) 
u 'i u 'j k 

and v^ v {r) = d 2 v/dr fl dr v (the indices fx and v denote space directions). Now we are 
back to a real problem of liquid theory, since we have only d degrees of freedom per 
molecule (the center of mass coordinates), and the number of clones, m, appears as a 
parameter in (fl3|). 

Once one has derived an expression for the replicated free energy, one can deduce 
from it the whole thermodynamics, as described above (Notice that the 'technical' 
approximation of neglecting the exchange of atoms between different molecules, as well 
as using a harmonic model, means that one really studies the IS in this computation, 
rather than the real free energy valleys). In all three cases, one finds an estimate of 
the Kauzmann temperature which is in reasonable agreement with simulations, with a 
jump in specific heat, from a liquid value at T > Tk to the Dulong-Petit value C = 3/2 
(we have included only positional degrees of freedom) below Tk- This is similar to the 
experimental result, where the glass specific heat jumps down to the crystal value when 
one decreases the temperature (Our approximations so far are similar to the Einstein 
approximation of independent vibrations of atoms, in which case the contribution of 
positional degrees of freedom to the crystal specific heat is C = 3/2). The parameter 
m*(T) and the cages sizes are nearly linear with temperature in the whole glass phase. 
This means, in particular, that the effective temperature T/m is always close to Tk, 
so in our theoretical computation we need only to evaluate the expectation values of 
observables in the liquid phase, at temperatures where the HNC approximation for the 
liquid still works quite well. 

A more detailed numerical check of these analytical predictions involves the 
measurement of the complexity, 

£ t = S(T) - S valle y(T) (15) 

The liquid entropy is estimated by a thermodynamic integration of the specific heat from 
the very dilute (ideal gas) limit. It turns out that in the deeply supercooled region the 



temperature dependence of the liquid entropy is well fitted by the law predicted in |48 
Sii q (T) = aT~ 2 / 5 + b, which presumably allows for a good extrapolation at temperatures 
T which cannot be simulated. As for the 'valley' entropy, it can be estimated as that of 
an harmonic solid. One needs however the vibration frequencies of the solid. These have 
been approximated by several methods, most of which are based on some evaluation of 
the Instantaneous Normal Modes (INM) [EH] in the liquid phase, and the assumption 
that the spectrum of frequencies does not depend much on temperature below Tk- 
Starting from a typical configuration of the liquid, one can look at the INM around it. 
In general there exist some negative eigenvalues (the liquid is not a local minimum of the 
energy) which one must take care of. Several methods have been tried: either one keeps 
only the positive eigenvalues, or one considers the absolute values of the eigenvalues 
16|, E7j. Alternatively one can also consider the INM around the nearest inherent 
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Figure 4. The configurational entropy versus temperature in binary mixtures of soft- 
spheres and of Lennard- Jones particles. The soft sphere result (left curve), from |45[ ], 
compares the analytical prediction obtained within the harmonic resummation scheme 
(full line), to simulation estimates of Su q — S va ii ey , where the valley entropy is that of a 
harmonic solid with INM eigenvalues projected onto positive eigenvalues (+), taken in 
absolute values (x), or taken around the nearest inherent structure (*). The squares 
correspond to the numerical estimate of the thermodynamic configurational entropy 
obtained by studying the system coupled to a reference configuration (see text, and Q 
for details). The Lennard- Jones result (right curve), shows as a full (black) curve the 
theoretical prediction obtained from the cloned molecular liquid approach jl6| ^|. The 
dotted (green) curve is the result from the simulations of [l7j and the dashed (red) 
curve is the result from the simulations of [fil) . Both simulations use the Su q — S va u ey 
estimate where the harmonic solid vibration modes are approximated by the ones of 
the nearest inherent structure. 



structure which has by definition a positive spectrum |45], |4(| [47], |51J . This procedure 
really measures the configurational entropy rather than the thermodynamic complexity. 
The computation of the thermodynamic complexity, using its definition as a system 
coupled to a reference thermalized configuration, has also been computed in |^5[ and 
turns out to be not very different from the configurational entropy, on the time and 
temperature scales which have been studied so far (they must differ on infinitely long 
time scales, as we discussed in the previous section). 

The results for the configurational entropy as a function of temperature are shown in 
fig.||, for binary mixtures of soft spheres and of Lennard- Jones particles. The agreement 
with the analytical result obtained from the replicated fluid system is rather satisfactory, 
considering the various approximations involved both in the analytical estimate and in 
the numerical ones. 



7. Conclusion 

Our knowledge on first principle computations of glasses is still rather primitive. 
Basically we have obtained, for the glass, the equivalent of the Einstein approximation 
for the crystal. Even within this simple scheme, doing the actual computation for the 
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glass turns out to be rather involved. What is most needed next is: on the analytical 
side, some better approximations of the molecular liquid state, allowing to go beyond 
the small cage expansions, and some reliable estimation of time scales in the regime 
T K < T < T c ; on the numerical side, some precise results in the glass phase at 
equilibrium fo^fl ; on the experimental side, some more measurements of the fluctuation 



dissipation ratio in the out of equilibrium dynamics. No doubt: "... there are still 
fascinating questions of principle about glasses and other amorphous phases..." f2j. 
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